Who is the Lord of the Rings: Majorana, Dirac or Lifshitz? 
The Spin-Orbit-Zeeman Saga in Ultra-cold Fermions. 
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We discuss the emergence of rings of zero-energy excitations in momentum space for superfluid 
phases of ultra-cold fermions when spin-orbit, Zeeman fields and interactions are varied. We show 
that phases containing rings of nodes possess non-trivial topological invariants, and that phase 
transitions between distinct topological phases belong to the Lifshitz class. Upon crossing phase 
boundaries, existing massless Dirac fermions in the gapless phase anihilate to produce bulk zero- 
mode Majorana fermions at phase boundaries and then become massive Dirac fermions in the 
gapped phase. We characterize these tunable topological phase transitions via several spectroscopic 
properties, including excitation spectrum, spectral function and momentum distribution. Since 
the emergence or disappearance of rings leads to topological transitions in momentum space, we 
conclude that Lifshitz is the lord of the rings. 
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Ultra-cold atoms have now become standard laborato- 
ries to test for existing or new theoretical ideas in atomic, 
condensed matter, nuclear and astrophysics. The major 
appeal found in these table-top experiments is the abil- 
ity to tune interactions, populations, species of atoms 
and dimensionality - which constitute the standard tool- 
box for investigations of interacting bosonic or fermionic 
systems. Very recently a new tool has been added to 
the toolbox: the ability to tune simultaneously spin- 
orbit and Zeeman fields in a system of ultra-cold bosonic 
atoms The same tool can also be used to study 
ultra-cold fermionic atoms [l|-[3| and to simulate differ- 
ent condensed matter systems such as topological insu- 
lators 0, non-centrosymmetric superconductors (Hf and 
non-equilibrium systems [sl , where spin-orbit coupling of 
the Rashba-type [7| is encountered. 

This direct connection to condensed matter physics 
inspired a new direction in ultra-cold fermionic atoms 
where spin-orbit coupling of the Rashba-type has been 
very recently investigated (8l-[l3|. However, spin-orbit 
fields currently realized in experiments involving ultra- 
cold atoms correspond to an equal superposition of 
Rashba 0] hi?(k) = vpi{—ky^^kxy) and Dresselhaus 14 1 
h2:)(k) = VD{ky:k + kxy) fields, leading to equal- Rashba- 
Dresselhaus (ERD) form [H, [l3| hERDi^) = vkxj, where 
vr = vd = v/2. Other forms of spin-orbit fields re- 
quire additional lasers and create further experimental 
difficulties [l5|, such that the current Zeeman-spin-orbit 
Hamiltonian created in the laboratory is 



Hzso(k) 



-hzCFz — hyCFy 



hERDiy)Cfy (1) 



for an atom with center-of-mass momentum k and spin 
basis It), U). The fields = -^r/2, hy = -S/2, and 
hERoO^) = vkx can be controlled independently and can 
be used to explore phase diagrams as achieved in ^^Rb 
experiments [1]. Here, Qr is the Raman coupling and S 
is the detuning. 



Hamiltonian: To investigate artificial spin-orbit and 
Zeeman fields in ultra-cold Fermi superfluids, we start 
from the Hamiltonian density H(r) = Ho(r) + l-Liir)^ 
with ?i = 1. The single-particle contribution is 



Ho(r) = ^tAj(r) f/f(k)l + Hzso(k)l ,^.'(r), (2) 

s,s' 



where K{k) = k^/(2m) — /i is the kinetic energy relative 
to the chemical potential /i, [H^^olkjjss/ are the matrix 
elements of the Zeeman-spin-orbit matrix Hamiltonian 
defined in Eq. ([T]), k = — iV is the momentum operator, 
and ^pl{r) creates fermions with spin s at position r. The 
interaction term 



ni{r) = -5V|(r)Vl(r)V;(r)^/'t(r), 



(3) 



is local and g represents a contact interaction strength. 
We define the total number of fermions as = TV-^ + N^, 
and the induced population imbalance as Pind = (^t ~ 
N^)/N. We choose our scales through the Fermi mo- 
mentum kE defined from the total density of fermions 
n = -\- = /c|./(37r^). This choice leads to the 
Fermi energy ce = kp/2m and to the Fermi velocity 
ve = kE/m, as energy and velocity scales respectively. 

We focus on the zero-detuning case 6 = {hy = 0), 
use the basis V^j(k)|0) = |k t), V^|(k)|0) = |k |), where 
|0) is the vacuum state, and write l-Lo{r) as the matrix 



Ho(k) = K{k)l - hzCTz - hERD{^)cry 



(4) 



The interaction Hamiltonian Hx(r) can also be converted 
into momentum space as H/(q) = — gh\c\)h{<c\) ^ where 
the pair creation operator with center of mass momentum 
q is 6t(q) = + q/2)V^j(-k + q/2), and g can 

be expressed in terms of the scattering length through 
V/g = -Vm/i4Tra,) + Ek l/(2ek). 

Helicity Basis: The matrix Ho(k) can be diagonalized 
in the helicity basis <l>|(k)|0) = |k ^), ^|(k)|0) = |k ^), 
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via a momentum-dependent SU(2) rotation. The helicity 
spins ^ and J| are aligned or antialigned with respect to 
the effective magnetic field heff(k) = h^z + hERDO^)^- 
The eigenvalues of the Hamiltonian matrix Ho(k) are 
e^(k) = K{k) - |he ff(k)| and ^^(k ) = K{k) + |heff(k)|, 
where |heff(k)| = y^h^ -\- /i|.^^(k) is the magnitude of 
the effective magnetic field. The interaction Hamiltonian 
T-Li{q) can be written in the helicity basis as 1-Li{q) = 

-9 J2a^^s ^1/3 (q) ^7^5 (q). where pairing is now described 
by the operator 



which has eigenvalues for the highest quasiparticle band 



BU^D = 5] A„;3(ki,k2)*l(ki)$t (k2) 



(5) 



and its Hermitian conjugate, with momentum indices 
ki = k + q/2 and k2 = — k + q/2. The matrix 
AQ,/3(k + q/2, — k + q/2) is directly related to the ma- 
trix elements of the momentum dependent SU(2) rota- 
tion into the helicity basis, and reveals that the center of 
mass momentum ki +k2 = q and the relative momentum 
ki — k2 = 2k are coupled and no longer independent. 

Tensor Order Parameter: From Eq. (j5j) it is clear that 
pairing between fermions of momenta ki and k2 can oc- 
cur within the same helicity band (intra-helicity pair- 
ing) or between two different helicity bands (inter-helicity 
pairing). For pairing at q = 0, the order parameter 
for superfluidity is the tensor AQ,^(k) = AoAc^^(k, — k), 
where Aq = —9^js{B^s{^))^ leading to components: 
A^^(k) = zAT(k)sgn [kx] for helicity projection A = +1; 
A^4(k) = -As{k), and A^^{k) = As{k), for helicity 
projection A = 0; and A^j|(k) = — iAT(k)sgn [/c^^] , for 
helicity projection A = — 1. The amplitudes AT(k) = 
Ao\hERDm/Mk)\ and As{k) = Ao/i./|heff(k)| re- 
fleet the triplet and singlet components of the order pa- 
rameter in the helicity basis. The Bloch-sphere relation 
|AT(k)p + |A5(k)p = |Aop, shows that the singlet and 
triplet channels in the helicity basis are not independent. 

Higher angular momentum pairing: In the triplet sec- 
tor A^^(k) and Aj|j|(k) contain not only p-wave, but 
also /-wave and higher odd partial waves, as seen from a 
multipole expansion of |heff(k)|~-^ = [h^ + h'^ERD^M)] 
for finite hz. Similarly in the singlet sector A^j|(k) and 
A^^(k) contain s-wave, d-wave and higher even par- 
tial waves, as long as the Zeeman field hz is non-zero. 
Higher angular momentum pairing occurs because the 
local (zero-ranged) interaction in the (t,i) spin basis is 
transformed into a finite-ranged anisotropic interaction 
in the helicity basis (^,^). 

Excitation Spectrum: The effective Hamiltonian in the 
helicity basis takes the matrix form 



Hsp(k) 









a(k) 



A;^(k) A;^(k) 

VA;,(k) A;,(k) 



AM(k) A^4(k) 









-a(k) J 



(6) 



El(k) = ^(a--^et + |A5(k)|2) +|AT(k)|2, 

and for the lowest-energy quasiparticle band. 



^2(k) = y(a_ + ^?2^ + |As(k)|2) +|Ar(k)p. 

The eigenvalues for quasihole bands are Es(k) = — £^2(k) 
and ^4(k) = -^i(k). The term a. = M^) - ^(k)] /2 
is the average energy difference between the helic- 
ity bands = -|heff(k)|; while the energy £,h+ = 
[^^(k) + ^^(k)] /2 is the averaged energy sum of the he- 
licity bands £^h+ = ^(k) = e\^ — fi. 

Notice that £;i(k) > £;2(k) > 0, but that only £:2(k) 
can have zeros (nodal regions) corresponding to the locus 
in momentum space satisfying the following conditions: 

a) £,h- = ~'\Jih+ 1^5(^)1^, which corresponds physi- 
cally to the equality between the effective magnetic field 
energy |heff(k)| and the excitation energy for the singlet 



component + [^^(k)^; and b) |AT(k)| = 0, corre- 

sponding to zeros of the triplet component of the order 
parameter in momentum space. 

Phase Diagram: Since only £^2(k) can have zeros, 
the low energy physics is dominated by this eigenvalue. 
In the ERD case, where \hERDO^)\ = v\kx\^ zeros of 
E2(k) can occur when kx = 0, leading to the follow- 
ing cases: (a) two possible lines (rings) of nodes at 
{k'^ + kl)/{2m) = /i + ^J}ll^^\A^ for the outer ring, 
and {kl + kl)/{2m) = fi - ^Jhl - |Aop for the inner 
ring, when /i^ — |Aop > 0; (b) doubly-degenerate line 
of nodes at (k^ + k1)/(2m) = /i for /i > 0, doubly- 
degenerate point nodes for ji = 0, or no- line of nodes 
for /i < 0, when /i^ — | Aop = 0; (c) no line of nodes when 
/i^ — |Aop < 0. In addition, case (a) can be refined into 
cases (a2), (al) and (aO). In case (a2), two rings indeed 
exist provided that ji > ^/h^^^\Ao\^. However, the in- 
ner ring disappears when /i = ^/h^^^JAo\^. In case (al), 
there is only one ring when < ^/h^^^\Ao\^. In case 
(aO), the outer ring disappears at /i = —\/h1 — |Aop, 
and for ji < —^/h^^^\Ao\^ no rings exist. 

In Fig. [TJ we show the phase diagrams of Zeeman field 
hz/^F versus interaction parameter l/{kFCis) for spin- 
orbit coupling v/vf = (a) and 0.28 (c), as well as 
induced population imbalance Pind versus l/{kFCis) for 
v/vf = (b) and 0.28 (d). We label the uniform su- 
perfluid phases with zero, one or two rings of nodes as 
US-0, US-1, and US-2, respectively. Non-uniform (NU) 
phases also emerge in regions where uniform phases are 
thermodynamically unstable. Possible NU phases in- 
clude phase separation, modulated superfluids and su- 
persolid. The US-2/US-1 phase boundary is determined 
by the condition ja = ^Jh\ — |Aop, when \hz\ > |Ao|; the 
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US-O/US-2 boundary is determined by the Clogston-like 
condition \hz\ = |Ao| when /i > 0, where the gapped US- 
phase disappears leading to the gapless US-2 phase; 
and the US-O/US-1 phase boundary is determined by 
/i = — |Aop, when \hz\ > |Ao|. Furthermore, 

within the US-0 boundaries, a crossover hue between an 
indirectly gapped and a directly gapped US-0 phase oc- 
curs at /i = 0. 
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FIG. 1: (color online) Phase diagrams of hz/ep and Find 
versus l/{kFCLs) for ERD coupling v/vf — (a), (b) and 
v/vf — 0.28 (c), (d). Uniform superfluid phases are labeled 
as US-0 (gapped, either directly or indirectly), US-1 (gapless 
with one ring of nodes), and US-2 (gapless with two-rings of 
nodes). The NU label describes the region where uniform 
superfluids are unstable. 

Dirac and Majorana fermions: Changes in nodal struc- 
tures of the order parameter are associated with bulk 
topological phase transitions of the Lifshitz class as noted 
for p-wave 3, and d-wave 18, 19| superfluids. Such 
transitions are possible here because spin-orbit and Zee- 
man fields induce higher angular momentum pairing in 
the helicity basis. In the US-1 and US-2 phases near the 
zeros of £^2(k), quasiparticles have linear dispersion and 
behave as Dirac fermions. The disappearance of nodal 
regions (rings) correspond to annihilation of Dirac quasi- 
particles with opposite momenta. The transition from 
phase US-2 to indirect-gap i-US-0 occurs through the 
merger of the two-rings at the phase boundary followed 
by the immediate opening of the indirect gap at finite 
momentum. However, the transition from phase US-2 to 
US-1 corresponds to the disappearance of the inner ring 
through the origin of momenta, and the transition from 
phase US-1 to the direct-gap d-US-O corresponds to the 
disappearance of the last ring also through the origin of 
momenta. 

The last two phase transitions are special because the 



zero-momentum quasiparticles at these phase boundaries 
correspond to Majorana zero energy modes if the phase 
(f{k) of the spin-orbit field hERD{k) = \hERD{k)\e''^^^\ 
where (fCk) = sgn [kx] 7r/2 and the phase ^(k) of the or- 
der parameter Aq = |Ao|e*^*^^^ hold the relation at zero 
momentum: (p{0) = —0{0) [mod(27r)]. This can be seen 
from an analysis of the eigenfunctions 

corresponding to the eigenvalue £^i(k) (with i = 2,3) 
where Uij^i^ = t^ij(k) are the elements of the unitary 
matrix U that diagonalizes the Hamiltonian Hgp. The 
emergence of zero-energy Majorana fermions requires the 
quasiparticle (quasihole) to be its own anti-quasiparticle 
(anti-quasihole): ^j(k) = ^^^^(k). This happens at zero 
momentum k = 0, where the amplitudes Uii{0) = /7*3(0) 
and Ui2{0) = [7^4 (0), leading to the conditions /i^ = 
/i^ - |Aop, and (f{0) = -6>(0) [mod(27r)], showing that 
Majorana fermions exist only at the US-l/US-0 and US- 
2/US-l phase boundaries. In Fig. [21 we show the Lifshitz 
transition from US-1 to US-0 phase, where nodal (mass- 
less) Dirac Fermions in the US-1 phase become bulk zero- 
mode Majorana fermions at the US-l/US-0 phase bound- 
ary, and then massive Dirac fermions in the US-0 phase. 

The commonality between bulk Majorana fermions 
(found here) and surface Majorana fermions (found in 
topological insulators or superconductors) is that both 
exist at boundaries: bulk Majorana zero-energy modes 
may exist at the phase boundaries between two topo- 
logically distinct superfluid phases, while surface Majo- 
rana zero-energy modes may exist at the spatial (sample) 
boundaries of a topologically non-trivial superconductor. 

Lifshitz Transition: The transition between different 
superfluid phases occurs without a change in symmetry of 
the order parameter tensor AQ,/3(k) in the helicity basis, 
and thus violates the symmetry-based Landau classifica- 
tion of phase transitions. However a finer classification 
based on topological charges can be made via the con- 
struction of topological invariants 0, [2o| . The number 
of rings £ corresponds to the topological charge associated 
with the surfaces of zero-energy quasiparticle excitations. 
Thus, for the US-0 phase £ = 0, while for the US-1 and 
US-2 phases, i = 1 and £ = 2, respectively. 

Spectral Function: An important measurable quan- 
tity is the single-particle spectral density [21| As{uj, k) = 
— (l/7r)ImGss(ia; = + k) for spin s =t,i, which can 
be extracted from the diagonal elements of the matrix 



G{iuj, k) = A{iuj, k)I + B{iuj, k) • a. 



(7) 



A{iuj,k) 

'^^ - \K. 



where the scalar function is 
KikMiLof - K\k) - lAop] - Mk)\'^[tu; 
K{k)]}/D{iuj^'k)^ and the vector function components 
are Bx,y{iuj,k) = {/i^,2^(k)[|/ieff(k)p + |Aop - {ioj + 
i^(k))^]}/D(icj, k), for the transverse and Bz{iuJ,'k) = 
/i,{[|Mk)|' - |Aop - {ico + K(k)f]}/D{iLo,k) for the 
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FIG. 2: (color online) Excitation spectra Ei(k) in the 
(kx — 0,ky,kz) plane illustrating the Lifshitz transition: the 
shrinkage of Dirac rings (US-1 phase) into Majorana zero- 
energy modes (US-l/US-0 phase boundary) and emergence 
of massive Dirac fermions (direct-gap d-US-0 phase). 



longitudinal parts. Here, D{iuj^'k) = nj=i(^^ ~ ^jO^)) 
and a is the vector Pauli matrix. In Fig. [3l we show 
^s(a;,k) in the plane of momenta ky-kz with kx = 
and uj = revealing the existence of rings of zero-energy 
excitations in the US-1 and US-2 phases. In the left 
panels the spectral densities for the US-1 phase are 
shown for spin ^ (a), where the ring is brighter than 
the ring for spin | (d). In the middle panels {b) and 
(e) ^s(a;,k) at the US-l/US-0 phase boundary is shown 
revealing the Majorana zero-energy mode. In the right 
panels (c) and (/) the spectral densities for the US-0 
phase vanish at a; = 0, since this phase is fully gapped. 




FIG. 3: (color online) The zero-energy spectral density 
As{oj = 0,kx = 0, ky, kz) at l/{kFCis) = 1-0 and v/vf = 0.28 
is shown in (a) and (d) for the US-1 phase with hz/ep — 
1.75, in (b) and (e) for the US-l/US-0 phase boundary with 
hz/ep — 1.59, and in (c) and (/) for the direct-gap d-US-0 
phase with hz/tp — 1-44. 

Momentum Distribution: A spectroscopic quantity 
that is routinely measured is the momentum distribution 



,(k) =T^[A(zc^,k)±5,(ic^,k)], 



(8) 



where the +(— ) sign corresponds to spin t (i)- Since 
ns(k) depends only on the energy spectrum Ej{\<i) and 
its derivatives, it is an even function of momentum k. In 
Fig. m we show ns(k) for l/ikFag) = 1.0 and v/vp = 
0.28 at the US-1 phase with hz/ep = 1.75 (top panels), 
at the US-l/US-0 phase boundary with hz/ep = 1.59 
(middle panels), at the US-0 phase with hz/ep = 1.44 
(lower panels). The left-most (right-most) panels show 
the momentum distribution for spin t (|) at /c^ = versus 
ky (solid-blue line) and versus kx (dotted-red line). It 
is very important to note the discontinuity of ns(k) at 
the location of the ring of zero-energy excitations in the 
US-1 phase (top panels), its change in behavior as bulk 
Majorana fermions emerge at k = (middle panels), 
and the transition to a direct-gap d-US-0 phase (bottom 
panels), where the ring of nodes has disappeared leading 
to smooth momentum distributions ns(k). 
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FIG. 4: (color online) Momentum distribution ns(k) for 
l/{kFCis) = 1.0 and v/vf = 0.28 at the US-1 phase with 
hz/cF = 1.75 (top panels), at the US-l/US-0 phase boundary 
with hz/cF = 1.59 (middle panels), at the US-0 phase with 
hz/cF = 1.44 (lower panels). The left-most (right-most) pan- 
els show the momentum distribution for spin t (i) at /c^ = 
versus ky (solid-blue line) and versus kx (dotted-red line). 

Summary: We showed that the presence of simultane- 
ous Zeeman and spin-orbit fields induces higher angular 
momentum pairing in the helicity basis, and identified 
topological phase transitions of the Lifshitz class via the 
existence of: rings of nodes in the excitation spectra, 
Dirac quasiparticles, bulk Majorana zero-energy modes, 
and topological charges. Lastly, we characterized dif- 
ferent topological phases via experimentally measurable 
quantities such as the spectral function and momentum 
distribution, and concluded that Lifshitz is the lord of 
the rings. 
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